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A Problem  in  Self-Heating  of  a Spherical  Body 

by 

S.  M.  Genensky 
National  Bureau  of  Standards 

Abstract 

An  analytic  steady-state  solution  is  developed  for  a spherical 
body  in  which  heat  is  generated  according  to  a first  order  uni- 
molecular  reaction  law,  and  lost  at  the  surface  in  accordance 
with  Newton's  law  of  cooling.  The  temperature  within  the  sphere 
depends  chiefly  on  the  radial  distance  from  the  center,  but  also 
on  the  ambient  temperature,  the  surface  heat  loss  coefficient, 
and  the  material  properties  which  are  assumed  constant. 

Introduction 

The  analysis  described  in  this  paper  was  pursued  in  connection 
with  an  investigation  into  the  self-heating  of  fibrous  materials 
carried  on  at  the  National  Bureau  of  Standards. 

Comparison  of  the  results  of  this  analysis  with  the  experi- 
mental results  indicated  that  the  assumption  of  a single  first 
order  unimolecular  reaction  was  an  over-simplification  for  the 
material  investigated  experimentally.  However,  the  analysis  serves 
to  indicate  the  relative  importance  of  the  various  parameters  enter- 
ing into  the  problem,  and  may  be  helpful  in  indicating  the  kind 
of  data  that  would  be  of  greatest  use  in  further  investigations 
in  which  the  assumptions  may  be  applicable. 

A homogeneous  sphere  of  radius  B is  generating  heat , under 
steady-state  conditions,  in  accordance  with  a first  order  uni- 
molecular  reaction  law  Ae“pp  and  losing  heat  from  its  surface 
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to  an  atmosphere  at  constant  temperature  according  to  Newton's 

law  of  cooling.  The  temperature  T(r)  of  the  sphere  is  independent 

of  time  and  depends  upon  the  radial  distance  r from  the  sphere's 

center.  Within  the  sphere  transfer  of  heat  by  means  other  than 

conduction  is  considered  negligible.  The  thermal  conductivity 

k,  gas  constant  R,  rate  of  reaction-frequency  product  A and 

activation  energy  E of  the  material  of  the  sphere  are  assumed  to 

dT 

be  known  constants  as  are  the  temperature  Tq  and  gradient  ^ = 0 
at  the  center  of  the  sphere.  The  problem  is  to  find  the  surface 
temperature  of  the  sphere  Tg  at  r = B and  the  heat  transfer 
coefficient  h between  the  surface  of  the  sphere  and  a 
surrounding  atmosphere. 


Analysis 


Mathematically  the  problem  becomes: 


(1) 


T = Tq  at  r = 0 


(2) 


(3) 


= -hZ  TB-Ta_/ 


(4) 


Let 


Then  equations  (l)-(4)  become: 


d^V  ^ 2 ^ 


1 

+ Ce  y — 0 


1) 


(5) 
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V = Vo  = g 


0 


(6) 


dr^  ^ 


at 


n 


= 0 


. ^ BhR  /-^  7 

E ^ Ta_/ 


(7) 

(8) 


It  is  further  assumed  that  (a)  the  temperature  gradient 
within  the  sphere  exists  and  is  continuous  throughout  the  inter- 
val 0^  1 and  (b)  that  the  absolute  value  of  the  difference 

between  1(7^)  and  Tq  , is  small  in  comparison  to  Tq. 

Under  these  restrictions  a technique  presented  by  Chambre^  proves 
useful  to  observe  that; 


1 


1 ^ ^ 

V Vq-^^V 


- ^ Zl/(i-f  )_7  - ^ 


V. 


I€ 


(9) 


i=0 


where  AV  = Vn  - V 


Now  since  series  coverges  and  for  a sufficiently 


1 ^ i 

large  positive  integer  N,  7—  (7  is  a very  good 

^ Vo  o 

approximation  to  y. 

Therefore  (5)  may  be  written  approximately  as 

lb?  + 2 dv  ^ 

d)^2  7 dry 

Now 


N 


42  + a 2 ce-%  ^ h: 


o 


= 0 


(0  < 


v^< 


1) 


N A ^ V • 

21  (7^)^  = Z1  (1-  “ = (N+1)  - BqV  + B2V2  - 

i=0  i-0  ° 

+ (-1)%VN 


(10) 


N 


n: 


% “ (n-k)'  ki 


n=0 


V, 


k 


k = 1,2,  : . . ,N 


where 
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and  in  the  special  case  V = V 


(N+l)-BiV  + B2V2-  ...  +(-l)  Bp|V  =1  N = 0,1,2, 


Therefore  (10)  may  be  written 

^ + Ce-fy  e*7Z  ZB1V-B2V2  + -(-l)»BMvh7  _ 


= 0 


Since  (11)  is  analytic  and  regular  except  at  = 0 and 

n 


further  since  0 at  = 0 the  solution  of  (11)  may  be 

written  in  the  form: 


V = V(/|)  - 


3-il 


iY]^  (o^y|^i) 

Recalling  (6),  equation  (12)  at  yields: 


Vo  - V(0)  = 


_ TqR  _ 


E 


3-0 


and  further  differentiating  (12)  once  and  considering  (7) 
it  is  found  that 


aq  = 0 

Differentiating  (12)  twice,  substituting  these  first 
two  derivatives  into  equation  (11)  and  rearranging  terms: 

Oca 

ZI(l2+i)a  -r.l-2  = _De+  - BaVS  + ...  -(-1)%?”  _7 

i=2  1 \ o 

N±l 

where  D = Ce“  Vq 


and  the  aq  are  given  by 


ai 


iz.2  q _ ]\T  _ 

1 dZDe-*-v7ZBlV-B2V2  . - (-D'^'BnV”/  7 

(i^+i ) (i-2 ) d/Y"2 

i=2,3,4... 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 
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Using  the  coefficients  obtained  from  (13),  (1^),  and 
(l6)  and  substituting  them  into  equation  (12),  V can  be 
computed  for  0^ 

In  particular  for  7|^=  1 equation  (12)  gives  the 
dimensionless  surface  temperature  V(l)  and  thus  recalling 
that  T (/|)  = pV (/^)  the  surface  temperature  T(l)  which  equals 
Tb  is  easily  found.  Differentiating  equation  (12)  once 
and  evaluating  this  derivative  at  7|  = 1,  h can  then  be 
found  by  solving  equation  (8),  for  all  the  other  factors 
involved  in  this  equation  are  either  given  or  have  been 
computed . 

The  first  six  non-zero  coefficients  have  been  evaluated 
and  are ; 


^o 

C _ 1 
a2  - - ^ e V 


o 


^4  - 120  Vq2 

8 C _ 2 -r  _ 1- 
"6  = Iw8'7o3  ^ 15  Vo  3-/  " 

c5  ,-122c£  122C  10  7 - - 

as  - 51840  Vo^  - 63Vq2  ' 21  Vq  3 -5® 

^ q1 -4032  c3  8428  i486  C 460  7 5 

®10  4435200  Vq5  ^ 405  Vo3  ■ 135  Vo^  27  Vo  ' 9 -^e-  Vq 

Observe  that  (i  ~ l,2,3,o..)  are  zero.  This 

follows  from  ap  = 0 and  that  spherical  symmetry  of  the  problem. 

Figure  1 is  a plot  of  V(l),  as  a function  of  C and  1/Vq , 
where  C was  allowed  to  vary  over  the  range  10“l<j^C^10l9 
and  1/Vq  took  on  selected  values  in  the  range  5^1/V^lOO. 
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The  points  were  obtained  using 


10 


(17) 


over  its  range  of  applicability  for  the  values  of  C and  I/Vq 
within  the  above  limits. 

In  the  foregoing,  it  has  been  assumed  that  the  temperature 
at  the  center  of  the  sphere  Tq  is  known,  and  also  that  the 
physical  condition  indicated  by  equation  (8)  is  satisfied. 

The  more  common  steady  state  problem  of  finding  the 
temperature  at  the  center  of  a sphere  which  generates  heat 
according  to  a first  order  unimolecular  reaction  law  and  loses 
heat  according  to  Newton's  law  of  cooling  may  be  solved  using 
the  analysis  described  above.  In  this  case  E,  R,  k.  A,  h,  T^, 

Tb  and  B are  assumed  known. 

If  the  sphere  being  considered  really  satisfies  the 

assumptions  of  the  problem  then  its  temperature  must  satisfy 

equation  (1).  Since  equation  (1)  is  of  the  second  order  only 

two  boundary  conditions  can  be  preassigned . Further  since  the 

temperature  at  the  center  of  the  sphere  must  remain  finite,  one 

dT 

of  these  preassigned  boundary  conditions  must  be  = 0 at  r - 0. 
This  leaves  but  one  free  boundary  condition.  However  at  the 
surface  of  the  sphere  both 


(18) 


and 


T = Tb  at  r = B 


(19) 
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must  be  satisfied;  and  since  only  one  of  these  two  boundary 

conditions  may  be  formally  assigned,  it  must  be  assumed  that  the 

other  is  consistent  with  the  problemo  Therefore  Figo  1,  which 

represents  equation  (17) ? may  be  used  to  solve  this  problemo 

Since  R,  E,  k,  A and  B are  known,  both  C and  V(l)  (recall  that 

Tb  = ^-*V(1))  can  be  computed*  Thus  using  Fig*  1,  1/Vq  may  be 
K 

found  by  interpolation  provided  that  the  point  (C,  V(l))  under 

consideration  lies  within  the  region  of  applicability  of  equation 

(17)  as  indicated  on  Page  6*  Since  the  center  temperature  of 

E 

the  sphere  is  T^ , that  is,  ^ problem  is  solvedo 

Here  no  use  was  made  of  condition  (l8)  and  as  mentioned 
earlier  this  condition  must  be  satisfied  independently,  if  the 
assumptions  made  are  fulfilled*  Thus  the  degree  of  consistency 
between  conditions  (l8)  and  (19)  serves  as  an  indication  of  the 
applicability  of  the  assumptions  made  in  this  analysis* 
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